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Abstract 

>. 

Equations are found for exact g-functions corresponding to integrable bulk and boundary 
flows between successive unitary c < 1 minimal conformal field theories in two dimensions, 
confirming and extending previous perturbative results. These equations are obtained via an 
embedding of the flows into a boundary version of Al. Zamolodchikov's staircase model. 
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1 Introduction 



The exact g-function [TJ[2] is a powerful tool for the study of integrable boundary flows, allowing 
the results of, for example, OS] to be extended to situations where the bulk is not critical. The 
initial proposals of [JJ were restricted to cases where the bulk theory possessed only massive 
excitations, and their scattering was diagonal both in the bulk and at the boundary. Recently, 
in [5], equations were introduced to describe the exact g- function for the simplest case were 
massless bulk degrees of freedom persist even in the far infrared, namely the flow between the 
tricritical and critical Ising models. As mentioned in [5], these 'massless' (^-function flow equations 
can be obtained from a consideration of the so-called staircase models, and as a result are 
naturally embedded in a much richer set of flows linking the boundary behaviours of all of the 
unitary c < 1 minimal models of conformal field theory. In this paper we will provide some 
more details of this larger pattern, and in the process propose equations to describe (^-function 
flows between all neighbouring pairs of unitary minimal models. For all of these cases beyond 
the first (which is the already-discussed flow between tricritical and critical Ising models), bulk 
scattering is non-diagonal. However, the g-function equations have a simple form which naturally 
generalises previously-seen structure^]. 

The staircase connection naturally leads us to equations which describe two-parameter fam- 
ilies of boundary perturbations, special cases of which match the one-parameter flows found 
perturbatively in [9] , and which therefore also match the results for fluctuating geometries found 
in [10]. We expect that our more-general sets of flows from superpositions of boundary states can 
be generalised yet further, to describe integrable bulk and boundary deformations of boundary 
conformal field theories with arbitrary numbers of boundary parameters. At the level of the exact 
^-function equations the generalisation is rather clear, and will be indicated below. However we 
will leave the detailed investigation of this point for further work, as the two-parameter situation 
is already quite involved. 

2 The staircase model, in the bulk and at the boundary 

The staircase model was originally introduced by Al.B. Zamolodchikov in [IT]; various generali- 
sations can be found in |12H15| . Its S-matrix encodes the diagonal scattering of a single massive 
particle of mass M, and can be obtained by the analytic continuation of the S-matrix of the sinh- 
Gordon model to those complex values of the coupling constant where real-analyticity holds. At 
the level of the Lagrangian the meaning of this continuation remains somewhat obscure, but as 
an S-matrix theory the model appears to make perfect sense, and leads to a consistent picture of 
finite-size effects described exactly by thermodynamic Bethe ansatz (TBA) equations which are 
derived in the standard way from the S-matrix. Trading the analytically-continued sinh-Gordon 
coupling for a real parameter 6q, the staircase S-matrix is 



and this leads to the following TBA equations for the pseudoenergy e(6) for the system on a 
circle of circumference R : 



1 As we were writing this paper, an alternative derivation of the diagonal p-function flow equations of [I, 2,5 , 
was presented, in [pj (see also [7]). It will be worthwhile, and appears in most respects to be straightforward, to 
generalise the approach of [pj to cover our new equations, but as we feel the staircase aspect is of independent 
interest we have decided to leave this point for the time being. We should also mention that boundary g-functions 
in the staircase model were previously discussed, amongst other things, in [8], but since that paper predated the 
exact g- function results of [TJ, its conclusions for situations with off-critical bulk were not correct. 
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where r = MR, and 



L(0) = ln(l + e -«W) , MO) = -^T9 lXim ■ 



(2.3) 



The ground state energy of the system is then given by E{R) = — ^ c e g- (r), where the effective 
central charge c e fj(r) is 

Ceff ( r ) = JL f rcosh(O)L(0)d0. (2.4) 



7T Z 



For later use we define 

sinh(f + SfS) 



(x)(e) 



sinh (| - ^) ' 



— sin(-7rx) /(2tt) 
cosh(6>) — cos(7rx) 



and also 



1 



'^) yW > ~ 2vrcosh(0) 
so that the staircase kernel 4>s(&) is 

For large values of 9 , this function is localised about 9 = ±9q, as illustrated in figure [TJ 



(2.5) 
(2.6) 
(2.7) 
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Figure [TJ The staircase kernel (f>s(6) for 9 = 20. 



The TBA system (|2.2p therefore couples the values of the pseudoenergy e(0) near to 6 with those 
near to 6 ±9o, and the behaviour of the effective central charge c e g (r) depends crucially on how 
many times the interval [0, 0q] fits into the range [— ln(^), ln(-)] , beyond which the value of the 
pseudoenergy is dominated by the driving term rcosh# in (|2.2p . irrespective of its coupling to 
values taken elsewhere. Referring the reader to |lipi3j for further explanation, the net result is 
that c e ff (r) develops a series of plateaux, or steps, which become more pronounced as 9$ — > oo. 
(The top curves on figures [3M3H below show c e fr(r) for 9q = 60, by which point the plateaux are 
already quite sharply defined.) Indexing the steps by an integer m = 3, 4 . . . , the (m— 2) th step 
is found for —(m—2)9o/2 <C ln(r) — (m— 3)6>o/2, and on this step, c e ff(r) — > c m as 9q — > oo, 
where 

Cm. — 1 
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m(m+l) 

is the central charge of the unitary minimal model A4 m . (More precisely, this holds in a double- 
scaling limit: pick f and 9q with —(m—2)9o/2 < ln(f) < — {m— 3)9$/ 2, and set r = f p and 
^0 = P^o with p > 0; then lim p _ i . 00 c g(r, 9q) = c m .) Furthermore, in the crossover region 
ln(r) ?» — (m— 3)0o/2 between the M m and A4 m _i plateaux, the staircase pseudoenergy e(0) 
tends uniformly on suitably-shifted intervals to the pseudoenergies ej(0) which solve the system 
of TBA equations introduced in [16J to describe the 0i3-induced flow from M. m to A4 m -i, a 
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flow which had previously been found perturbatively in [171118]. (This interpolating theory was 
denoted by AdA^ in [16J.) This is illustrated in figure [2] for m = 4. 
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Figure [2] The function L{9) at various values of r, for 9q = 60. From the upper to the lower 
curve, ln(r) is equal to —50, —40, —30, —20 and —10. This range covers the crossover between 
tricritical and critical Ising models, and the form of L(0) around 6 = +30 and —30 approximates 
the functions ln(l + e _£l( - e_30 ^) and ln(l + e _e2 ^ +30 ^) from the corresponding interpolating TBA 
of US] (cf. figure 1 of 0). 



The natural interpretation of these results - also backed by perturbative studies, in the spirit 
of [17] . at large m [H] - is that there is a one-parameter family of integrable quantum field 
theories with renormalisation group trajectories which, in the limit 0q — ^ °°; approach the union 
of the renormalisation group trajectories of the MAm^ theories. From now on we will assume 
this to be the case, and use it to deduce equations for the flow of the (/-function in these same 
theories. 

To treat the boundary staircase model using exact ^-function techniques, we need a conjecture 
for its boundary reflection factor R(0). It is natural to suppose that this can be obtained through 
the same analytic continuation of the sinh-Gordon boundary reflection factor as yielded the 
bulk S-matrix (|2.1|) . The boundary sinh-Gordon model with no additional boundary degrees of 
freedom has a two-parameter family of integrable boundary conditions [20] . and its reflection 
factor follows from that of the first sine-Gordon breather, found in [21] (see also, for example, 
[22]). Further continuing to the staircase values of the coupling, this leads to the following 
reflection factor: 

_ i8p \ ( 3 , iOg \ 
V2/V4 27T A3 27T > 



R(0) 



(2 2X2 2X2 2X2 2) 



(2.9) 



where E and F are two parameters whose relationship with the original two parameters of the 
boundary sinh-Gordon model will not be relevant below. In the sinh-Gordon model E and F are 
often real, but for the staircase model it will be more interesting to consider them at the complex 
values for which real-analyticity is preserved, as with the continuation of the bulk coupling. 
Hence we set 

E = i?K F = ^ (2 . 10) 

7T 7T 

with Obi and Q\a real and, without loss of generality, non-negative. There is an obvious extension 
of this ansatz to incorporate n boundary parameters Obi ■ ■ ■ Obn '■ 



R(0) 



(IX! 



2tt A 



3 , iOa- 

4 2tt ■ 



-r-rn (I _ i6 bk w 1 . i8 bk 
llk=l\2 27T A 2 "r" 2?r 



(2.11) 



For n > 2 this does not correspond to an integrable sinh-Gordon boundary condition of the 
simple form treated in [20], but it can be realised by the addition of a stack of n— 2 defects next 
to such a boundary [23][24"]. 

For all of its subtleties at intermediate scales, the boundary staircase model in the far infrared 
is simply a massive diagonal scattering theory, both in the bulk and at the boundary. Its exact 
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g-function should therefore be given by the formula proposed in [TJ. Explicitly, for a boundary 
at the end of a cylinder of circumference r = MR, 



In g(r)= In g (r) + In g b (r) (2.12) 



where 



^o(r) = tij un i^k) ' ' ' TT%)^ + e ^ - ^ ■ ■ ■ - 0l) (2 - 13) 



n=l 

and 



]ng b (r) = \ [ d9 (0 6 (0) - 0s(20) - \ 6(6)) L(0) (2.14) 

z JR 

where e(6) solves the TBA equation L(6) = ln(l + e~<^), 5 (0) = -^W lnS ( 9 ) = 

0(0 — 6q) + 0(0 + 0q), and, restricting attention to the n = 2 case of (|2.1ip for simplicity, 

mo) = -~*m 

+ 0(9 - 9 bl ) + 0(9 + 6 bl ) + 0(0 - 9 b2 ) + 0(0 + fe2 ) . (2.15) 

As in [5] , this normalisation for <j) b differs by a factor of two from that used in [TJ[2] ; we also used 
0(0) = —0^1^(0). These equations are straightforward to implement numerically, and can also 
be treated analytically in various limits. In the next section we report some of the results of 
this analysis for the full staircase model, while in section H] we show how the staircase g- function 
equations decompose in suitable limits into sets of equations which govern exact ^-function flows 
in the interpolating theories AiAm^ ■ 

3 The boundary staircase flows 

Following the split (I2.12D of In g into the sum In g + In g b , the single- integral piece In g b naturally 
splits into three further terms as In g b = In g bl + In g b2 + In g bs , with 

lng bl = l [ d0(-0(0)-I<5(0))L(0), (3.1) 

z JR 

lng b2 = lJ R d9 (0 ( 3 } (0 - ±0 O ) + ( | } (0 + ±0 O ) - 0(20 - O ) - 0(20 + 9 ))l(9) , (3.2) 
lng b3 = \ ! dO (0(0 - 6 bl ) + 0(0 + 6 bl ) + 0(0 - 9 b2 ) + 0(0 + 62 ))L(0) . (3.3) 

1 JR 

Notice that In g bl and In g b2 are independent of the boundary parameters, and mg bl only depends 
on the bulk parameter 9q implicitly, via the function L(9). 

In the large-0o limit, the full g- function passes through a series of plateaux as lnr varies, 
its value on each plateau always matching a (conformal) (7-function value, or a product of such 
values, for the conformal field theory seen by the bulk theory at that value of r. Some features 
of this behaviour can be seen in figure [3l where for simplicity the values of 9 b i were chosen such 
that their associated boundary transitions always coincide with bulk transitions. These and 
other aspects will be analysed in more detail later. 



4 




In(r) ln(r) 
Figure EJ:: 9 bl = 120. Figure EJi: 9 bl = 180. 

Figure [3) Staircase g- function flows for various values of 6 b \. In each plot 9q = 60, and the value of 
9 b 2 is scanned from to 200 in steps of 2. To aid the eye five ^-function flows have been highlighted. 
The lowest-lying of them has 9 b 2 — 200; reading from left to right, this is then joined by the flows for 
6 b 2 = 150, 100, 50 and respectively. The short horizontal lines (light blue online) show the logarithms 
of the conformal g-function values (I3.5[) for the 'pure' (Cardy) boundaries of each minimal model visited 
by the staircase flow. The other g-function plateaux correspond to superpositions of these boundaries, 
as discussed in the main text. 
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To give some precise formulae, we start by recalling some facts about the minimal model 
M. m and its conformal boundary conditions. The model has central charge c m = 1 — m /^ +1 ) ; 
and the permitted irreducible highest-weight representations of the Virasoro algebra are labelled 
by pairs of integers a and b with 1 < a < m— 1 and 1 < b < m subject to the identifications 
(a, b) ~ {m— a, m+1— b). They have highest weight 

{{m+l)a - mb) 2 - 1 
Kb = 4m(m+l) ' (3 ' 4) 

In particular the bulk field (pis , which induces the interpolating flow from M m to .M m _i , has 
scaling dimension 2/113 = 2(m— l)/(m+l). The most general conformal boundary condition is a 
superposition of 'pure' (Cardy) boundary conditions [25]. There is one such boundary condition 
for each irreducible highest-weight representation (|3.4p . and its boundary entropy (or g- function) 

ism 



\ sin(^)sin( bT 



l) 



g(m,a,b) = [— — r , (3.5) 

Notice that this formula respects the identification (a, b) = (to— a, m+1— 6), and has the ad- 
ditional symmetries g(m,a,b) = g(m,m—a,b) and g(m,a,b) = g(m,a,m+l—b). This is the 
Z 2 'spin flip' symmetry of minimal model boundary conditions which in the case of the Ising 
model maps between spin up and spin down Dirichlet boundaries. For superpositions of Cardy 
boundaries, the boundary entropies simply add. 

We will also need some more detailed information about the form of the function L(9), 
illustrated for some sample values of ln(r) in figure [2j Suppose the bulk theory is near to the 
minimal model A4 m , so that ln(r) satisfies 

- (m-2)e /2 < m(r) < -(m-3)0 o /2 • (3.6) 

Setting a = 21n(l/r) — (m—3)9o, we have < a < flo and, starting from 6 w ln(l/r), L{9) 
exhibits an alternating series of plateaux of lengths a and 6q — a, the i th plateau being centred 
at 8 = Zi, where 

Zi = (m-2-*)0 o / 2 J i = 1,2, ...2m-5. (3.7) 

(For later comparison with the TBA systems for AiAm^ it will be convenient to count these 
plateaux starting from the right.) The plateau values of L{6) can be found as explained in [13] . 
Adapting slightly the notation from [16, 26J , define constants x a and y a by 



sin 2 !^) sin 2 (^ 

V m+l / \ m , 

= — 4 f, l + Va= 4 ( • (3.8) 

sin^^l sin 2 ^ 



Starting from 9 = +00, L(9) is close to until 9 ~ ln(l/r), and then has height ln(l + x%) for 
ln(l/r) - a < 9 < ln(l/r), then ln(l + y 2 ) for ln(l/r) - 9 <C 9 < ln(l/r) - a, then ln(l + x 3 ) 
for ln(l/r) — 9q — a 9 <C ln(l/r) — 9o, and so on, before returning to for 6 <C — ln(l/r). In 
full, the plateau values of L(9) are 

\n(l+x a ) : z 2a -3 - «/2 < 6» < z 2a -3 + «/2 , a = 2...m-l; (3.9) 

ln(l+y a ): z 2a _ 2 - {0 -a)/2 « < z 2 a-2 + (9 -a)/2 , a = 2 . . . m-2 , (3.10) 

and the complete sequence between 9 = — ln(l/r) and 9 = +ln(l/r) is 

{ In(l+x m _i), ln(l+y m _ 2 ), ln(l+x m _ 2 ) , ln(l+y m „ 3 ) . . . ln(l+y 2 ), ln(l+x 2 ) } . (3.11) 
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As a shorthand we will refer to the intervals in the set (|3.9p as x-type, and those in the set (|3.10p 
as y-type. From (|3.8p . ln(l+xi) = ln(l+yi) = 0, and so we can formally add these two constants 
to the end of the sequence (|3.1ip while remaining consistent with the values taken by L(9) in 
the corresponding intervals, and likewise add ln(l+y m _i) and then ln(l+x m ) to the beginning. 
The symmetries x a = x m+ i_ a , y a = y m -a reflect the more general symmetry L(9) = L(—8). 

With these preliminaries completed we return to the exact y-function g(r). Three parts of 
lng(r) do not depend on the boundary parameters: lng , lng bl and lny 62 . These functions only 
undergo transitions at the values of ln(r) where there is a bulk crossover, that is at ln(r) = 
— (to— 3)#o/2, Tn = 3, 4 ... . The effective equations governing these transitions in the large-#o 
limit will be treated in the next section; here instead we will suppose that r satisfies (]3.6p so 
that the bulk theory is close to the minimal model Aim- Then g , g bl and g b2 are approximately 
constant, and given by the following formulae: 



In g (r) 



In 



In 



m(m + 1} 



K m(m + 1) 



sin 



(m— 1)tt 
2m 



sin • 



sin 



m+l 



2(m+l) 

sin — sin 

m 



m+l 



for m odd 



for to even 



(3.12) 



(3.13) 



ln 9w (r) = -U(0)= < 



|ln(l+x 



(m+l)/2j 



In 



iln(l+?/ m /2 



■I In 



sin 



sin 



2 TT 



m+l 



for to odd 



for m even 



(3.14) 



(3.15) 



ln ft> 2 ( r ) 



|ln(l+y ( 



m-l)/2 



■Jin 



sin 



2 {m-Vj-K 
2m 



sin 



2 TT 



for to odd 



(3.16) 



±ln(l+x m/2 



These results are exact in the limit {9q — > oo, 



■iln 



TO- 



8111 



2(m+l) 



sin 



for to even . (3.17) 



m+l 



-2)0 o /2 < ln(r) < 



TO- 



-3)6» /2}. The 



formula for hig will be derived in section [4] below. Those for lny fel and lng fe2 follow from the 
fact that (f)(9) and <f>( x }(9) are only significantly non-zero near to 9 = 0. This means that the 
integrals (|3.ip and (|3.2p only receive contributions from the regions 9 ~ and ~ #o/2, where 
L(9) is, by (13. 9p and (I3.10p . approximately constant. Pulling L(9) out of each integral and using 



(W 



d0< 



d9 0(3)(0) 



1 

4 



(3.18) 



together with the plateau values given by (|3.9[) and (|3.10p leads to (|3.14p - (]3.17p . 

When lny , lny fel and ln<7 b2 are summed, the pieces which depend on whether to is odd or 
even cancel, leaving the following simple result, valid for all values of to: 



lng +lng bl +lng b2 = In 



TO(m + 1, 



7T 7T 

sin — sin 

TO TO + 1 



(3.19) 



This is the logarithm of 5(771,1,1) or g(m,m— 1,1), the boundary entropy of the conformal 
boundary condition associated with the bulk vacuum field or its Z2 spin flip conjugate. In fact 
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it is not surprising that this partial sum should be equal to the logarithm of the full boundary 
entropy for some boundary condition, as will become clear as we examine the behaviour of the 
remaining part of In g, namely ln<7 fe3 . 

Since In g b3 depends on the boundary parameters, we would expect it to undergo transitions 
not just where the bulk crossovers occur, but also, possibly, at energy scales related to pure- 
boundary transitions, and this turns out to be the case. The integral (|3.3|) receives contributions 
from 8 w ±6 b i and 8 ps ±8 b2 ; if these regions lie within the x- and y-type intervals (|3.9p and 
(I3.10P then, pulling L{8) outside the integrals and recalling that L(9) = L(—6), 

\ng b3 {r) = \{L(8 bl ) + L{8 b2 )) (3.20) 

and the value of In g b3 (r) will not change for small changes in r. Conversely, In g b3 {r) will undergo 
a crossover whenever r is such that either 8 b \ or 8 b2 lies on a boundary between the intervals (13. 9ft 
and ()3.10p . This means that there will be boundary transitions associated with the parameter 
8 b \ at 

ln(r) = -k9 - 8 bl , A: = 0,1,... (3.21) 

and 

ln(r) = -ke + 8 bl , k = A,A+1, . . . (3.22) 

where A = \28 b \ / 9q] , the smallest integer greater than or equal to 26*^1 / 80 • An analogous formula 
holds for 8 b2 . 

If ln(l/r) is smaller than 9 b \, L{8) is effectively zero near to 8 = ±8 b \ and the term ^L(8 b \) 
ceases to contribute to In g b3 (r). This explains why the plots in figure [3] stabilise with increasing 
8 b x, with the parts of the plots with ln(r) > —9 b \ being independent of b \. If ln(l/r) is smaller 
than 8 b \ and 8 b 2, then In g b3 (r) is zero and the logarithm of the full g-function is given by the sum 
In <7 (r) + \ng bl (r) + ln g b2 (r), which we already observed was equal to In g(m, 1, 1) when the bulk 
theory is on the plateau corresponding to M. m - If the limit 8 b \ — > 00, 8 b2 — > 00 is taken before r is 
varied, we see that there is one flow which simply moves through the (1, 1) boundary conditions in 
the successive minimal models, its g- function being given by In g{r ) = In g Q (r ) +ln g bl (r ) +ln g b2 (r ) 
for all values of r. (This is why the partial sum (|3.19p is itself the logarithm of a boundary 
entropy.) For ln(r) > —180, this flow is matched by the lowest-lying curve of figure [3H. Notice 
that since \n.g b3 {r) is manifestly positive, all other g- function flows must lie above this limiting 
curve, an off-critical generalisation of the fact that at a fixed point the lowest-possible boundary 
entropy is always found for the (1, 1) boundary condition. 

For smaller values of 8 b \ and 8 b2 the picture becomes more complicated, as can already be 
seen from figure O Nevertheless it is still possible to formulate general rules for the boundary 
conditions which are visited. The plateau behaviour of L(8) means that the sequence of boundary 
conditions seen for any given 8 b \ and 8 b2 depends on the intervals (|3.9p and (|3.10p that they (and 
their negatives) find themselves in as ln(r) varies. From (|3.1ip and the immediately-following 
remarks, the possible values of L{8) on these intervals are the elements of the set 

{ln(l + ^),ln(l+y b )} (3.23) 

where the indices a and b lie in the ranges 

a G {l,--- ,m}, b G {1, • • • ,m- 1} (3.24) 

and the value of 0, found for \8\ 3> ln(l/r), arises when a is equal to 1 or m, or b is equal to 1 
or m—1. The symmetries x m +\- a = x a , y m ~b = Ub could have been used to restrict the indices 
a and b to 

in ~\~ 1 txx — 1 

a e {!,■■■ , — - — }, 6e{l, — - — } for m odd 

m 

and a, b G {1, • • • , — } for m even (3.25) 
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but for reasons to be explained below it will be convenient to keep with the larger ranges. A 
given pair of boundary parameters corresponds to two (possibly equal) values of L(9) from the 
set (|3.23p . and, via the exponential of (|3.20p and (|3.19p . to a value for the g- function. We found 
that this value can always be expressed as a sum of Cardy g- function values (|3.5p . according to 
the following rules, where we introduce the convenient notation [x a ,y b ] and so on to denote the 
particular combinations of boundary conditions which arise: 

L(6bi) L(9b2) Boundary condition 

ln(l+x a ) ]n(l+y b ) [x a , y b ] = (b, a) ^ ^ 

ln(l+x p ) ]n(l+x q ) [x p ,Xg] = (1, \p-q\ + l)k(l,\p-q\+3)k ■ ■ ■ k(l,m-\p+q-m-l\) 

ln(l+y r ) In(l+y a ) [y r ,y s ] = (|r-a|+l, l)&(|r-a|+3, 1)& • • • k(m-l-\r+s-m\, 1) 

with the same result for 9 b \ o 9 b2 . Notice that if all indices are restricted to the reduced ranges 
(|3.25p . the rules simplify to 

L(6 b \) L{0 b 2) Boundary condition 

ln(l + x ) ln(l + y b ) (b, a) (3 27) 

ln(l + x p ) ln(l + x q ) (l,|p-g|+l)&(l,|p- 9 |+3)&---&(l,iH-g-l) 

ln(l + y r ) m(l + y s ) (|r-s|+l, l)&(|r- s|+3, 1)& • • ■ &(r+s-l, 1) . 

The identifications of (/-function values with specific boundary conditions implied by (13.26P 
should be treated with care for a couple of reasons. First, ambiguity arises from the equalities 
g(m,a,b) = g(m,m—a,b), g(m,a,b) = g(m,a,m+l—b). These are related to the symmetries 
L(9) = L(-9), x a — x m _|_i_ a , y a — Vrn—a °f L(9) and its plateau values, and, more precisely, to 
the following symmetries of the boundary condition combinations given in f|3.26|) : 



[x a ,yb] = [x m +i-a,y m -b] [x a ,yb] = [x m +l-a,yb] = [Xa,y m -b] 



[Xp,Xq\ — [x m _|_l_p, 2^m+l— q] [Xpj Xq\ — [^m+1- pi %q\ — [Xpi Xm+1 — q] (3.28) 



ym-s 



[yr,y s ] = [y m - r ,y m - s ] [y r ,y s ] = [y m ~r,y s ] = [y r ,Vr, 

where the overbar denotes the Z2 'spin flip' symmetry mentioned just after (|3.5p . acting on 



individual Cardy boundaries as (a, b) = (m—a,b) = (a,m+l— b) . As remarked in [5], it should 
be possible to resolve such ambiguities in a systematic fashion by studying the renormalisation 
group flows of the inner products of finite- volume excited states [27)128] with the boundary states, 
but we shall leave this for future work. Second, there are sometimes further, more accidental, 
degeneracies in the set of non-negative-integer sums of Cardy g-function values - for example, 
g(5, 1,3) = 2g(5, 1, 1), so that in A4$ the (1,3) and (1, 1)&(1, 1) boundary conditions cannot be 
distinguished by their ^-function values alone. Nevertheless, and modulo the spin flip ambiguity 
just described, (|3.26p is the only set of decompositions we have found which works in a uniform 
fashion for all m. From now on we shall assume that it is correct, and mostly leave the spin flip 
ambiguity implicit. 

Continuing to suppose that the bulk theory is in the vicinity of the bulk fixed point M. m , we 
now let ln(r) vary from the lower to the upper end of the range (j3.6l) . that is from — (m— 2)9q/2 to 
— (m— 3)9q/2 . The centres of the intervals (|3.9p and (|3.10p remain fixed, at 9 = Zi, % = 1 . . . 2m— 5, 
but the widths of the x-type intervals, a = 21n(l/r) — (m— 3)9q, decrease from #0 to zero, while 
those of the y-type intervals, #0 — «> increase from zero to 9q. Thus so long as \6 b \\ < (m— 2)9^/2 
and 9 b \ is not an integer multiple of 9q/2, the regions 9 ~ ±#&i move from x-type intervals to 
y-type intervals during this process, and the value of L(±9 b i), and hence that of the g- function, 
undergoes a change. For brevity we will phrase the rest of the discussion in terms of the intervals 
seen by L(9 b i), but we could equally look at L{— 9 b i). By the Z2 ambiguity just discussed this 
might lead to different boundary conditions being assigned but since the ^-functions are blind to 
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their difference, we will ignore this issue for now. In fact, it will be convenient to allow 9b\ (and 
6*62) to take positive and negative values, so all options are in any case covered. 

Suppose, then, that at the start of the process 9bi is in the x-type interval centred at 9 = Z2 r -3, 
and #52 is hi the x-type interval centred at 6 = Z2 S - 3 , corresponding to plateau values for L{9b\) 
and L(9b2) equal to ln(l + x r ) and ln(l + x s ) respectively. The y-type interval that Obi moves 
to depends on its position relative to 22r-3- If &bi > -22r-3 then L(9b\) moves to the plateau 
ln(l + y r _i), while if Obi < Z2r-3 it moves to the plateau ln(l+y r ). An identical set of possibilities 
occurs for L(9b2), with a transition which may occur before of after that in L{0b\) depending on 
the relative distances of 9bi and 9b2 from the centres of their original (x-type) intervals. Putting 
these ingredients together gives the following 'skeleton' of transitions from an initial situation 
where L(0bi) = ln(l+x r ) and L(9b2) = ln(l+x s ) : 

[yr,y s ] < [x r ,y s ] > [y r -i,y s ] 



[y r ,x s ] i [x r ,x s ] > [y r _i,x s ] 



[yr,y s -i] < [x r ,y a -i] > [y r -i,y s -i] 

This diagram encapsulates the boundary flows seen while the bulk is in the vicinity of M m , if 9bi 
and 9b2 are in the intervals [z2 r -3 — 0o/2,Z2 r -3 + 0o/2], [z2 S -3 — Oq/2, Z2 S -3 + #o/2]. The complete 
set of flows seen near to M m is the union of a number of such diagrams, so as to cover the full 
ranges \0bi\ < (m— 2)0q/2, |^ 2 | < (m— 2)9q/2. To give one example, at m = 5 this results in a 
7x7 grid, the lower (9b2 > 0) half of which is: 



[y4,x 3 ] < — [x 4 ,x 3 ] — > [y3,x 3 ] < — [x 3 ,x 3 ] — > [y2,x 3 ] < — [x 2 ,x 3 ] — > [yi,x 3 ] 



[2/4,2/2] < [^4,2/2] > [2/3,2/2] < [X3,y2\ > [2/2,2/2] < [X2,V2] > [2/1,2/2] 



[2/4,^2] < [X4,X 2 ] > [2/3,^2] < [x 3 ,X 2 ] > [2/2,X2] < [x 2 ,X 2 ] > [yi,X 2 } 



[2/4,2/1] < — [£4,2/1] — ► [2/3,2/1] < — [%3,yi] — > [2/2,2/1] < — [x2,yi] — > [2/1,2/1] 

The entries in square brackets can be converted into specific boundary conditions by using 
the dictionary (|3.26p . to give the picture shown in figure U] below. We note once more that these 
predictions are made modulo the Z2 ambiguity in the relationship between ^-function values and 
boundary conditions. The options chosen here, which follow from the rule formulated earlier, 
are consistent with predictions made in, for example, 0E9], but we have not attempted to 
confirm them using exact g-function techniques. As already mentioned, this would require the 
computation of inner products of states other than the ground state with the boundary state. 

Notice that the figure is symmetrical about the diagonals 9b\ = 9b2 and 9b\ = —9b2, while 
negating either Obi or 9b2 individually has the same effect as the Z2 spin flip. This second feature 
means that the boundary conditions in the middle column on the figure, 9bi = (or equivalently, 
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the top row shown, 9 b 2 = 0) are mapped into themselves under the spin flip. As follows from 
(|3.28p . the corresponding properties hold for all other values of m, if it is assumed that the rules 
(|3.26p are correct in general. 
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Figure 31 A diagram of the pure boundary flows appearing near to the bulk fixed point VW5, for 6q = 60 
and various values of 6 b i and 9^2- The slightly-thicker vertical and horizontal arrows (red online) indicate 
the flows and boundary conditions which are mapped to themselves under the Z2 spin flip. 



To complete the picture we must examine what happens when the bulk theory moves between 
two neighbouring fixed points, say M m and M m _i. The bulk flow occurs when ln(r) varies 
through a region centred on ln(r) = — (m— 3)9q/2 which is of size of order one as 9q ~~ ^ °°- As 
this happens the form of L{9) changes: the x-type plateaux for A4 m shrink to zero size, while 
the y-type plateau of height ln(l + y a )\M m , a = 1 • • • m— 1, for M. m becomes the x-type plateau 
ln(l + x a )\M m ~i f° r M-m-l, with y-type plateaux for M. m -i opening up between these as ln(r) 
increases further. (By (|3.8p . ln(l + y a )\M m = l n (l + Xa)\M m ~n so there is no sudden change in 
heights implied by the redesignation of y-type to x-type plateaux over the transition.) 

Now consider the behaviour of the logarithm of the full g- function, Ing = lng + hig fol + 
ln<7 b2 + hiy 63 , as M. m flows to M. m ~\- From (|3.19p . the sum of the first three terms, lny + 
ln<? 61 + lny b2 , changes from In g(m, 1, 1) to In g(m— 1, 1, 1). The behaviour of the remaining 
piece, hiy b3 , depends on the values of 9 b \ and 9 b2 . Suppose first that both \9 b i\ < (m— 3)#o/2 
and \9 b 2\ < (m— 3)^o/2, so that neither boundary parameter has become decoupled at the point 
of the bulk transition under discussion. There are then three cases: 

(i) If neither 9 b \ nor 9 b 2 are at the centre of what had been an x-type plateau for M m , that is if 
0&1,0&2 7^ (m-3)0 o /2- k9 , fe = 0,l,...m-3 (3.29) 

then both 9 b \ and 9 b 2 will have undergone a 'pure-boundary' transition of the sort described 
earlier, from an x-type onto a y-type plateau, before the M. m — > M. m -\ transition is reached. 
Thus the conformal boundary condition seen just before the bulk transition occurs will correspond 
to some pair [y r , y s \\Mm where 2 < r, s < m— 2, corresponding to 9 b \ and 9 b 2 lying in the intervals 
((m-3)9 /2-(r-l)9 , (m-3)0 o /2-(V-2)0 o ) and ((m-3)9 /2-(s-l)9 , (m-3)9 /2-(s-2)9 ) . 
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After the transition the plateau values will not have changed but their interpretations will have, 
to the pair [x r , x a ]\Mm-i • Translated into specific conformal boundary conditions using (|3.26p 
at m and m — 1 the flow is therefore 

(/,1)&(/+2,1)& ••• k(g,l) M m 

I (3.30) 
(l,/)&(l,/+2)& ••• k(l,g) M m -l 

where 

/ = \r— s\+l , g = m— 1— \r+s— m\ and 2<r,s<m—2. (3.31) 

Via (|3.28p and the symmetry under 6{,i o 9^2 the full set of options is explored by restricting r 
and s to the fundamental domain 2<r<s<m— 2, r + s < m. In fact, (|3.31|) is equivalent to 
/ and g in ()3.30p being restricted by 

1 </<<?< m-1, f-g€2Z. (3.32) 

Notice that these, the generic flows, always start at 'sinks' on networks of pure-boundary flows 
such as figure HI and end on 'sources' on the corresponding network one minimal model down. 

(ii) If one of Obi or 9^2 lies at the centre of an x-type plateau for M m , then it remains on that 
plateau right up to the moment of the bulk transition, after which it will instead lie in the centre 
of a y-type plateau for M m -\. If this centre is located at 9 = (m—3)8o/2—(s—2)8o, 2 < s < m— 1, 
then the corresponding value of L(9) moves from ln(l+x s )|_M m to ln(l+y s _i)|^t m _ 1 . The other 
plateau value simply changes its designation from \n(l+y r )\M m t° hi(l+x r )|_M m _ 1 , as in case (i). 
Thus the boundary condition flow is [x s ,y r ]|x m — > [y s -i,%r]M m -ii or 

(r, s) M m 

I (3.33) 
(s-l,r) M m -i 

where 

2<r<m-2, 2 < s < m-1 . (3.34) 

(iii) Lastly, if both O^i and 9^2 are at the centres of x-type plateaux for M. m , say at (m— 3)9q/2 — 
(r—2)9o and (m — 2>)9q/2 — (s—2)9q with 2 < r, s < m—1, then reasoning as above the boundary 
condition flow is [x r ,x s ]\ Mm -> [y r -i, y s -i]\M m -i, or 

(l,/)&(l,/+2)& ••• &(1 )5 ) M m 

I (3.35) 
(/, 1) & (/+2, 1) & ■■■ & (ff-2,1) Mm-i 



where this time 



or equivalently 



/ = |? — s|+l , g = m— \r+s— m— lj and 2 < r, s < m— 1, (3.36) 



l<f<g<m, f-g€2Z. (3.37) 



For these (least-generic) cases the flows are always from sources to sinks on neighbouring pairs of 
pure-boundary networks such as figure HI and decrease by one the number of superposed Cardy 
states. 

Finally we must treat the cases where either one or both of \9bi\ and is larger than 
(m— 3)9q/2. Then the corresponding plateau values of L(9) simply flow from zero to zero. If 
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in addition neither Obi nor 0^ lie at the centre of an x-type plateau, then it is easily seen 
that the situation is covered by case (i) above, if the indices r and s are allowed to take the 
additional values of 1 and m—1 (recall that ln(l+yi)|x OT = ln(l+y m -i)|x m = ln(l+xi)|x m _ 1 = 
ln(l+a; m _i)|^f m _ 1 = 0). Thus the combined story is that there are flows of the form (|3.30p for 
every pair (/, g) satisfying 

l<f <9<m-l, f-ge2Z, (3.38) 

(Taking f = g = l or f = g = m-l gives the flows (1, l)|x m -> (1, l)\M m -i and (m-l, l)\M m -> 
(l,m—i)\M m -t which occur when both \0bi\ and \0b2\ are larger than (m—3)0o/2, so that lng b3 
remains zero throughout the flow.) Last of all, if, say, \0b2\ > (m— 3)#o/2 while Obi is at the 
centre of an x-type plateau, then the flow is \yi,x s ]\M m -> [xi,y s -i}\M m -i or [Vm-i , %s] \M m ~> 
[x m -i,ys-i]\M m -i with 2 < s < m—1, which simply means that the indices in ()3.33p can be 
given the enlarged range 

\<r<m-l, 2<s<m-\. (3.39) 
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Figure[5j A depiction of the flows corresponding to figure|3ji, for which B\>\ — 180 and = . . . 200. 
The rows are labelled by the corresponding bulk minimal model; nodes are then labelled by the 
boundary condition within that model. The highlighted flows, labelled in bold along the top of the 
figure, correspond to the values of Qbi that were also highlighted in figure [Sh 1 . Although some flow 
lines appear to cross, this is an artifact of the projection onto the page, and does not occur in the 
full multidimensional space of flows. 



These rules can be combined to understand the sequences of g- function flows seen in figure [3] 
and further illustrated in figures [5] and [6l Consider Obi = 180, 0b2 = 50, one of the highlighted 
flows in figures[3H and[5j Focussing on the part of the flow beginning at M$, the initial g-function 
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value is close to that of the boundary condition (1,4). Since 9f,i > ln(l/r) everywhere in this 
part of the flow, only 8 b 2 has an effect on the subsequent trajectory. The centre of the L{8) 
plateau associated with the boundary condition (1,4) at Ms is at 8 = 30, and since 8^2 > 30 
the flow within A4$ is to (3, 1). Then as the bulk theory flows from A4$ to Mj, the boundary 
condition flows to (1,3). The centre of the associated L{8) plateau is then at 8 = 60, and since 
8^2 < 60, the flow within .M7 is then towards (3, 1), and so on. Repeating this exercise for other 
values of 8 b 2 leads to the set of flows illustrated in figure the second highlighted flow of which 
corresponds to 8 b 2 = 50. Note that the pure-boundary flows within A4§ on figure [5] match the 
flows on the bottom row of figure HI 



50 100 150 200 




Figure[6] A depiction of the flows corresponding to figurc[3b, for which Obi = 60 and 0b2 = . . ■ 200. 
Other labelling is as on figure [5] 

We remarked while discussing figure [3] that the g- function plots stabilise as 8 b i increases. 
The same feature can be seen on comparing figure [5] with the equivalent diagram for 8 b \ = 60, 
figure [6j For .A/f 3 and .M 4 the boundary conditions and flows appearing in figure [6] match those 
in figure El but for the higher minimal models the boundary conditions appearing are different. 
This corresponds to how the plots in figures [3b and [3jd coincide for In r > —60. 

4 Exact (/-functions for MA { + ] 

Finally, we are ready to obtain the effective equations which govern the ^-function flows for the 
boundary versions of the bulk interpolating theories MAm^ ■ These theories can be thought of 
as bulk perturbations of the minimal models A4 m by their ^13 operators, with the sign of the 
perturbation chosen so that the infrared limit of the model is A4 m -i, the next minimal model 
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down. Our predictions for the associated boundary flows can be read from the results of the last 
section, but to obtain exact equations, we need to take a careful limit of the g-function formula 
in parallel with the double-scaling limit of the bulk TBA equations. 

To fix notations we first review the situation for the bulk TBA. The TBA system proposed 
in [16] for M.Am ^ involves m — 2 pseudoenergies e\ . . . e m _2, coupled together by the following 



system of TBA equations: 

ei{ e) = he e - [ <t>{e-e')L 2 {e')de' 

e a (9) = - f </>(e- 9')(L a ^(9') + L a+1 (9')) d6' a = 2... m-3 
Jr 

e m - 2 (0) = \re~ e - [ 0(9 - 9')L m ^(9') d9' (4.1) 

where L a (9) = ln(l + e~ €a ^) , (f)(9) = l/(27r cosh(0)) is as in (|2.6p . and r sets the crossover scale. 
The effective central charge is then 

ceff(r) = (e e L x (9) + eT 6 L m ^(9)) d9 (4.2) 

and as ln(f) increases through a region of size of order 1 about the origin, c e fi(r) moves from c m 
to c m _i, consistent with the corresponding bulk flow being from A4 m to M. m -\. 

These same equations emerge from a double-scaling limit of the staircase model if we set 

ln(r) = -(m-3)0 o /2 + hi(f) (4.3) 

and then take the limit 9q — > oo keeping ln(r) finite, before finally allowing ln(f) to vary from — oo 
to +oo. In particular, the pseudoenerg ies € a (9) are recovered from the staircase pseudoenergy 
e(9) in the 9q ^ oo limit by setting 

e a (9) = e(9+ (m-l-2a)9 /2) , a = l...m-2 (4.4) 

and only allowing 9 to vary over the full real line after the limit has been taken. 

It turns out that in this limit the staircase (^-function formulae can be rewritten in terms 
of the limiting pseudoenergies e a (9) and various constants which can be calculated in terms of 
the plateau values of the staircase pseudoenergies. We start with the infinite series part of the 
g- function, f|2. 13|) . Crucial to the analysis is the double-bump shape of the kernel 4>s(&), shown in 
figure [H which causes each multiple integral contributing to the sum to localise onto a collection 
of subregions of W 1 . In each of these subregions, the staircase pseudoenergy is either constant, 
or else is uniformly well-approximated by one of the interpolating- flow pseudoenergies e a (9). 
Rewriting the formula for In g in terms of these constants and functions leads to the effective 
equations which govern the (/-function flow in A4Am^ ■ 

More precisely, since 4>s(9) = 4>(6—6o) + <ft(9+8o), each term </>s , (#i + 6>2)</>s'(#2 — #3) ■ ■ ■ 4>s(Sn — 
9\) in the sum in f|2. 13j) can be expanded as sum of 2 n terms of the form 

0(0i + 9 2 - ai0 o )0(#2 - 3 - a 2 9 ) ■ ■ ■ <\>(9 n -9 X - a n 9 ) (4.5) 

where each = ±1. The decay properties of (f)(9) mean that it is only non-zero for 9 ~ 0, so 
for the above term to be non-zero as #0 ~> 00 we require 

01 + 6» 2 w ai 9 

02 — 63 ~ «2#0 



a n 9o (4.6) 
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from which it follows that Bu « Qu, k = 1 . . . n, where 



V 



h \ 

h 



' n 
h J 



/! 1 



2 00 



1 1 
1 1 
-1 1 



\ / "1 \ 

«2 



V i 



-1 -1 



-1 



(4.7) 



1 / V «n / 



and the integral over R n has indeed localised, to a set of 2 n regions of size of order one as 9$ — > oo, 
which become infinitely separated in this limit. The coordinates of the centres of these regions 
are either all even multiples of Qq/2, or all odd multiples of 0o/2, depending on whether n is even 
or odd. 

For each region of integration we must consider the behaviours of the 'measure factors' 
1/(1 + e e ^ fc ^), k = 1 ... n. As for the function L(9) discussed above, these factors exhibit a series 
of plateaux interleaved by transition regions, at 8 ~ 9 a = (m—l—2a)0o/2, a = 1 . . . m—2. Within 
these transition regions, the measure factors are well-approximated in the 9q oo limit by the 
functions 1/(1 + e €a ), by (|4.4p . In between these regions the measure factors are approximately 
constant, and can be expressed in terms of the numbers y a \M m = x a\M m -n a = 1---W— 1. 
Taking these considerations into account, the terms in the sum in (I2.13P fall into two categories: 

(a) If m + n is odd, every 8k satisfying \9k\ < (m— 3)9/2 lies in a transition region for e(9), so 
that its measure factor remains nontrivial even after the 9o oo limit has been taken. We 
denote the part of In go consisting of these these terms by In g A (f) . 

(b) If m + n is even, every 9^ lies inside a plateau of e{9) after the 9q — > oo limit has been taken, 
so that the corresponding measure factor becomes constant. We denote the (f-independent) 
part of In go consisting of these these terms by In g B . 

For the (a) terms, the values of e(9 ~ 9k) vary as f varies, and only reach plateau values in 
the UV and IR limits, these values being Xk+i\M m m the UV and Vk\M m -i m t ne IR. Rewriting 
the formulae in each subregion of integration in terms of the limiting pseudoenergies e a (9) using 
(|4.4p and shifting the integration variables to remove all appearances of 9o, we found that In g A 
can be rewritten in the 9q — > oo, r finite limit as 



\ng A (r) = V / antiT^n^A^)^, 



n-V 



l )cf>(9 n +9 l ) (4.8) 



m+n odd 



where the (m — 2) x (m — 2) matrix A{9) is given by 
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1+enW 

l 

l+e E 3(9) 
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l +e £ m-3W 





/ 

(4-9) 

and antiTr (K), the anti-trace of an M x M matrix K, is defined as the sum of its anti-diagonal 
elements, or equivalently 



antiTr (if) = Tt(KJ) , where Jjj = S^m+i- 



(4.10) 
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The measure factors for the (b) terms are by contrast constant throughout the relevant 
integration subregions. They can therefore be pulled outside their integrals, leaving only the 
various factors of (f)(0) . This leads to the following expression for \ng B : 

\ng B = V -UntiTr(5 n ) f 0(0i-0 2 )0(0 a -03) • • • 0(0„~i-0„)0(0„+0i) dB v ■ ■ d6 n 

»>1 iU jR " 

m+n even 

= E J^^i^) = E Tr (B n J) (4.11) 

n>l n>l 
m+n even m + n even 

where the tridiagonal (m — 3) x (m — 3) matrix I? is equal to the limit as Oq °° of 

/ o — r ^- T ••• 

1 o . ••• u 



rh=z \ ••• 



1+e 



o o ••• — , x m _ 6 , o — T .. 



m — 6 £ 



V ••• r4i=i7T 

The explicit form of this matrix can be found using the L(9) plateau values (|3,lip and is 







l-^rirw- 

sin 2 — 

o ••• o o o 



V 



••• 1 





sin 2 — _ sin 2 — 

777 || 1 771 

■ o (m-3)ir U ^ . 2 (m-3)7T 

sin i — sin j — 



or, more concisely, 



/{l+Xa+l) (4.12) 



where x a = x a \j^[ m _ 1 = S1 ^i^j^ — 1, and l a b is the incidence matrix of the j4 m _3 Dynkin 
diagram. Note that B is the transpose of a matrix which arises in the analysis of small fluctuations 
about stationary solutions of an associated Y-system [16] and has eigenvalues 

A fc = 2cos(— ), k = 2,3...m— 2. (4.13) 

m 

For later use we note that [B, J] =0, and furthermore that the eigenvector ipk of B corresponding 
to the eigenvalue satisfies 

JV fc = {-l) k ^k ■ (4.14) 
As explained in appendix A, this information is enough to evaluate (14. IIP in closed form, 
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with the result 



ln 5z?l A f r 




- — - 2m for m odd (4.15) 

V 3 ^ 




'MA},. 

/i \ 
/ / o \ z 1 \ 

for m even . (4-16) 

A missing piece of the staircase discussion from the previous section can now be filled in, namely 
the formulae f|3. 12|) and f|3. 13j) for the value of \n.g Q \j^i m when the bulk staircase theory is in the 
vicinity of M. m - In terms of the limiting ^-function equations for A4Am \ this number is equal 
to the UV limit of In g A (r) + lng B as f — > 0. Considering the limiting forms of the matrix A(6), 
given by (|4.9j) . as f — > 0, it is straightforwardly seen that 

lng A (0)\ MA w = lng B \ MA ^ (4.17) 

and so 

^g \ Mm = ]ng B \ MA (+^+]ng B \ MA M (4.18) 

which, via (jOB"]) and (l4"46l) . leads immediately to ff37T2]) and (13431) . 

The remainder of the exact g-function, In g b = In g bl + In g b2 + In g fe3 , is more straightforward 
to analyse. In the limit {0o °°; hir finite}, the behaviours of the first two terms, lng fol and 
ln<7 fe2 , depend on whether m is even or odd. 

The presence of 4>(9) and 5(9) in the expression ()3. 1|) for lng bl means that this term is 
determined by the staircase pseudoenergy close to 9 = 0. When m is odd this remains nontrivial 
as the double-scaling limit is taken and, using the pseudoenergies defined in (14. 4p . In ^r bl becomes 

]ng bl = -i / d9 (<f>(0) + ±<5(0))ln(l + e"^ W ) , (4.19) 

JR 

which has UV and IR limits given by (13. 14H at 7W m and (13.15P at Ai m -i respectively. When m 
is even, the staircase pseudoenergy instead becomes constant near to 9 = 0, giving 

ln 06i = ~\ ln (! + VtuMmJ ( 4 - 2 °) 

matching f)3. 15j) at M. m and also (|3. 14p at M. m -\. 

In contrast, the formula (|3.2[) for mg b2 involves 0(3) (0 ± ^0 O ) and 0(20 ± 0o), meaning that 
lng fe2 is determined by the behaviour of L(9) close to 9 = ±0o/2. For m odd, the limiting form 
of L(9 « ±0o/2) is constant, and 

ln<7 fe2 = -|ln(l + i/2^i|M m ) (4.21) 

as in (|3.16p at M. m and (|3.17p at M. m -i- When m is even, L(9 ±0o/2) remains nontrivial and 
in the limit 

ln 5fe2 = / d0(<W0) - 0(26»))ln(l + e"^ (e) ) , (4.22) 

with UV and IR limits given by ()3.17p at M. m and ()3.16p at M. m -\ respectively. 

To allow the last part of ln<7 fe , \ng b ^ to retain a non-trivial r-dependence in the limit, we 
pick two integers a\ and 02 with < a, < m— 1, write the boundary parameters 0fti and 6'f,2 as 

bi = 1(771-1-2^)00 + (4-23) 



18 



for i = 1,2, and then take the 9q —> oo limit keeping 6 b \ and b2 finite. Given the specification 
(|4.4p of the effective pseudoenergies e a (8), for 1 < Oj < m— 2 the staircase expression (|3.3p for 
ln<7 fe3 then reduces to 



2 /■ 

i=i ^ 
2 r 

= V / ^^-^)ln(l + e -^ (9) ) (4.24) 
i=i ® 

where the symmetry e a (#) = e m _i_ a (— 0) of the ground-state pseudoenergies for .MA^ was used 
in going from the first line to the second. If either a, is equal to or m— 1, then the staircase 
pseudoenergy diverges in the region of 9 bi , and the corresponding term in (|4.24|) is zero in the 
MA& limit. 

To summarize the results of this section, our final expressions for the two-parameter families 
of exact g-functions for A4Am\ indexed by a pair of integers a± and a 2 and expressed in terms 
of the rescaled variables f, b \ and 9 b2 , are as follows: 

• m odd: 



lng aia2 (f,§ bl ,e b2 ) =l n ^iy^nZj - \ jf d0 (0(0) + i<5(0))ln(l + e ^ W ) 

+ ln# A (r) + ln 5b3 aia2 (f, M , b2 ) 5 (4.25) 



m even: 



^g aia2 (f,9 bl ,6 b2 ) = ln^^y s in^j + jf M^^G) - 0(20)) ln(l + e 

+ ln 5A (f) + ln 5fe3 aia2 (f, 4l, 4 2 ) • (4.26) 

In both cases 5^4 is given by (14. 8p . the sum over n running through even integers for m odd, and 

odd integers for m even, with the pseudoenergies involved solving the bulk MAm ' TBA system 
(|4.ip . The term g b3 is as defined in (I4.24p . The constant terms result from adding ln^g, 
given by (|4,15p or (|4.16p . to lng b2 for m odd, and to mg w for m even. The remaining integral 
term is In g bl for m odd and In g b2 for m even. Formally setting one or both of a\ and a 2 equal to 
or m—1, as discussed after (|4.24p . incorporates the limiting one- and zero- parameter families 
of flows found by deleting the 9 b \ and/or 9 b2 dependent parts of hig b3 . If both are deleted so 
that ln(7 fe3 is identically zero, then by the results of the last section the (/-function flow should be 
from g(m, 1,1) at r = to g(m—l, 1, 1) as r — > 00. This can be checked directly. Simplest is the 
UV limit. For m odd, using f|3. 14|) and (|4.17p . the first three terms on the RHS of (|4,25p tend to 




m 



sin —tt / \ v 7 a /sin 



1 



m+l I / \ y m+1 



which is equal to g(m, 1,1). Similarly, for m even, (|3.17p and (|4.17f) imply that the first three 
terms on the RHS of (j4.26j) become 





^ -1 in 1 ^"T iiHi^i r"'~: ,j 1 <«^2s, 



and the value of g(m, 1,1) is again reproduced. 
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In the f — > oo, IR, limit the matrix A(9i) defined in (|4.9p again becomes independent of 
in the central region ~ where the integrals in (|4.8p have their support, and tends to the 
following (m— 2) x (m— 2) matrix: 







sin 









sin 





2 2tt 
m— 1 


















\ 






In terms of Air , 






ln# A (oo) 








£ 

n>l 

m + n odd 








2 (m — 3)71 



2 n+2 



Tr((A m ) n J) 



/ 

(4.29) 
(4.30) 



To evaluate this sum, observe that the central (m— 4) x (m— 4) sub-matrix of Air is the matrix 
B for .M-A^j]*^. Furthermore, the zero values of all entries in the first and final rows means that 
the central (m— 4) x (m— 4) sub-matrix of (.Air)™ is just (I?| 



(+) 



) n , and the entries on the 

first and final rows are again zero. Hence the traces and antitraces of (Ani) n are equal to those 
of (B\ M . (+) r,and 



m#(u(°°) 



'MA 



(+) 



In 50 



s \MA 



(+) 



Mi 



(+) 



allowing (|4.15p and (|4.16p at m— 1 to be used to compute lngt^ 00 ) 
(|3.15p and (|4.16p at m— 1, the first three terms of ([4.25P therefore become 

1 



(4.31) 

. For m odd, using 



In 



4V 

— Wsm — 

m / V ?7i 



In 



sm 



+ In 



m 



1 



m— 1 



sm 



TO— 1 



which is indeed equal to g(m—l, 1, 1). Similarly, for m even, in the r — > 00 limit the first three 
terms of (14261) are, by pTBll and (14151) at m-1, 



In 



2 \ 4 / vr 
— Wsm — 
m / V 



I In 



sm 



(m-2W 
2(m-l) _ 



sm 



+ In 



TO— 1 



m 



1 



I sin / (!Z^K 



2(m-l) 



(4.32) 



sm 



TO— 1 



which again matches g(m—l, 1, 1). 

More generally, the full equations (|4.25p or (|4.26p predict a collection of two-parameter fam- 
ilies of flows, indexed by the two integers a\ and 02- The relevant calculations have been carried 
out in the last section and we won't repeat them here. Instead, in figure [7] below we show a 
typical family of flows, where r = a\ + 1 and s = 02 + 1. All flows for Obi and #&2 finite start 
in the far UV at the [x r ,a; a ] boundary of M. m , a superposition of Cardy boundaries given by 
the rule (|3.26p . If both O^i and 0^2 are zero, the flow is directly downwards to the [y r —i,ys-i] 
boundary of M. m -i, driven by the bulk perturbation. Nonzero values of Obi and 0b 2 correspond 
to the addition of boundary perturbations to the bulk perturbation, and cause the trajectory to 
visit other boundaries on its way from UV to IR, as can be read from the figure. If either of a\ 
or 02 is equal to 1 or m—2, so that one or both of r and s is equal to 2 or m—1, the cube shown 
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in figure [7] truncates, the equalities x\ = y\ or £ m _i = y m -2 within M. m -\ causing one or two 
rows on the bottom face of the cube to fuse together, resulting in the flow patterns illustrated 
in figure The same phenomenon was seen earlier for the staircase model, and is reflected in 
the diagonal lines of flows running from top right to bottom left in figures [5] and 



A 




Figure [71 A cube of AiAm^ flows, where r = a.i+1 and s = ci2+l. Flows on the outer faces occur for 
(9;,i — > ±oo and/or #f,2 — > ±oo; those on the top face are within the minimal model A4 m , and those on the 
bottom are within A4 m _i. The central vertical axis corresponds to Obi = ®bi — 0. Conformal boundary 
conditions at fixed points of the flows are labelled according to the scheme given in equation (|3.26j) . 




Figure El Truncations of the cube of flows for r = 2 (left) and r = s = 2 (right), corresponding to 
setting either a\ = 1 or a\ — a,2 — 1 in (|4.24l) . Other details are as on figure [Jj 
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5 Conclusions 

In this paper we have shown how exact methods can be used to study interpolating boundary 
flows in two-dimensional integrable models in situations where both bulk and boundary are away 
from criticality. Our results for flows between minimal models have confirmed and extended 
previous perturbative studies. In addition we have shown how these flows can be embedded into 
larger manifolds of boundary integrability via the staircase model, generalising the picture seen 
in the bulk. In fact the staircase description is remarkably economical - the bulk S-matrix fj2. If) 
and the boundary reflection factor (12. 9p (or (12. lip ) together encode not only all of the unitary 
c < 1 minimal models, but also all of their Cardy boundary conditions and a variety of their 
superpositions, once fed into the general TBA and exact (^-function machinery. 

There are many directions for future work. Some perturbative checks of the exact equations 
for the tricritical to critical Ising bulk and boundary flows were undertaken in [5], but further 
tests higher up the series would be valuable, as would a more detailed study of the two-parameter 
boundary flows discussed in section and their extensions via defects to incorporate further pa- 
rameters through the more general reflection factor (|2.1ip . The formulae found in section U] 
describe, for the first time, exact off-critical ^-functions in situations where the underlying scat- 
tering theory is non-diagonal, and it will be interesting to generalise the approach of Pozsgay [6] 
to cover such cases. At the same time, many other multiparameter families of integrable models 
with nontrivial intermediate scaling behaviours are now known, including generalisations of the 
staircase models [T2HT5] . and the Homogeneous Sine-Gordon (HSG) models [30H32] , Hence there 
is plenty of scope to obtain more elaborate exact ^-function flows using the approach adopted 
in this paper. Finally, it is noteworthy how the embedding of non-diagonal bulk and boundary 
scattering theories within higher-dimensional manifolds of integrability achieved by the staircase 
and HSG models manages to 'abelianise' their TBA descriptions. It would be very interesting 
to know how general this phenomenon is, and to understand it at a deeper level. At the very 
least, it demonstrates once again that simple exact S-matrices and reflection factors can hide a 
great deal of internal structure. 
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Appendix A Proofs of equations (14.151) and (14.161) 

To derive (I4.15P and (|4.16p from (|4.1ip . we first convert the final sum, over traces of products 
of powers of B with J, into a sum of traces of pure powers of related matrices, after which the 
identity 



will allow the sum to be evaluated in terms of the eigenvalues of these related matrices. There 
are two cases. 





n=l 
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• If m is odd, then the sum in (|4.11f) is over odd powers of B. Shuffling indices and using the 
symmetries of B it can be checked that B n J = (BJ) n for all odd n, and so 



hi 



9b = ~^£~(§^) 



n 

n>l 

nodcl 



_ ll„ Det ( 1 + p J ). (A.2) 
8 Det(l-iBJ) v ; 

From (|4.14p . the eigenvalues of BJ are (— l) k Xk, which for m odd are the numbers 2 cos(27r/m), 
2cos(47r/m) . . . 2 cos((m— 3)7r/m), all with multiplicity two. Hence 



Det (1 + 


\BJ) 


A TV 4 

cos — COS 

rn 


2tt . 
m 


4 (m— 3)-7r 
• • COS 2m 


Det (l - 


\BJ) 


■ 4 7T '4 

sir ^_ sin* 


2?r , 

m 


■ 4 (m— 3)7r 
' • 8111 2m 



(A.3) 

To evaluate this we call upon some well-known trigonometric identities. For the numerator we 
use 

7r 2tt (n - l)ir 1 
cos — cos — • • • cos = —, — -T77T for n odd (A. 4) 

n n 2n 2( n " 1 )/ 2 V ; 

and for the denominator 

i 2f (n — l)ir n 
sin — sin — • • • sin = T (A. 5) 

which implies 



sin — sin — • • • sin = for n odd. (A. 6) 

„ „ o„ on— 1 v ' 



7T 2 27T 2\ n ~ I' 71 " n 

— sin — • • • sin 

n n 2n 

Using (|AT4l) and (|X6]l . 

— \ 2 

2m I 



Det (1 + \BJ) 1 sin 4 _ / 4 sin 4 ^i^' 



(A.7) 



Det(l-iSj) m 2 cos 4 (™_dlZL \m sin 2 ^ y 
and so (|A,2p is in agreement with (|4.15p . 

• If m is even, then the sum in (|4.1ip is instead over even powers of B. This time we make use 
of the identity, valid for m even, that 

Ti(B n J) = Tr(B n - B' n ) (A.8) 

where B' is the (m — 3) x (m — 3) matrix which only differs from B in the following elements: 

B m -2 m-4 = B ' m -2 m = 0. (A-9) 

2 ' 2 2 ' 2 

To prove this identity we write the LHS as 

Tr(B n J) = B ili2 Bi 2 i 3 ■ ■ ■ B in ^i n B inm ^i x -2 (A. 10) 

and note the following properties of B: 

B iti _i = B i i+ i, all other entries 0; (A. 11) 

Bij = B m _2-i,m-2~j ■ (A. 12) 

The first of these means that (|A. 10|) can be interpreted as a weighted sum over all n-step paths 
on the ^4m-3 Dynkin diagram which start and finish at pairs of conjugate nodes i\ and m—2—ii, 
i\ = 1 . . . m— 3, and move by one link at each step. Likewise Tr(B n ) and Tr(B' n ) are weighted 
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sums over n-step paths on the same Dynkin diagram, but which this time start and finish at the 
same node where again i\ is summed from 1 to m—3. 

These observations imply that both sides of (|A.10p are zero for n odd, and so from now on 
we can take n to be even (which is the case of direct interest in the current context). Due again 
to (jA.lip . a term in Tr(B n J) with i\ < must include the element B m-2 m , and a term with 

i\ > must include the element B m-2 m-4 . Assume first that i\ < and consider 

Z 2 2 Z 

Bi^Bi^ ■ ■ ■ Bi n _ 1 i n Bi nj7n -2-ii (A. 13) 

for some particular %i . . . i n . Suppose B m-2 m appears for the final time at Bii , which means 

2 2 p p+ 

that Bi p+1 i p+2 = Bm m+2 . By (|A.12j) . the value of (|A. 13|) is unchanged if all indices i q with q > p 
are replaced by their conjugates m — 2 — i q . The resulting term appears in the trace not of B n J, 
but of B n . Similarly, every term in the expansion of Tr(B n J) with i\ > can be equated 
with a term in the expansion of Tr(i? n ) with i\ > . Finally, the terms in the expansions 
of Tr(S n J) and Tr(i? n ) with i\ = ^j^- are already equal, since i± is then equal to m — 2 — i%. 
Thus Tr(B n J) is equal to the sum of the terms in the trace of B n that include either B m-2 m 

2 2 

or B m-2 m-4 at least once. Now the trace of B' n as defined by (1A.9I) is equal to the trace of B n 

2 2 

minus the terms where B or B m — 2 77?. — 4 appears at least once. Hence Tr(B n — B' n ) gives 

2 2 2 2 

the required terms and (|A.8P holds. 
If m is even we therefore have 

^9b = E ^T2 T ^ Bnj ) 



n even 



\ E ^(m n -(W)") 



_ (A , 4) 

using (jA.lj) . To find the eigenvalues of B' ', Det(B' — XI) can be expanded about the middle row to 
see that the characteristic polynomial of B' is proportional to A times the product of the (equal) 
characteristic polynomials of the upper-left and lower-right (m— 4)/2 x (m— 4)/2 submatrices of 
B', which we denote B[ and B' 2 . These submatrices can be fully diagonalised using the (m— 4)/2 
eigenvectors ipk of B with k odd: from (|4.14|) . these eigenvectors satisfy Jipk = — ipk- Hence their 
middle components are zero, while the neighbouring two are the negatives of each other. Given 
the definition (|A,9p of B' this means that projecting each ip^ for k odd onto its first (m— 4)/2 
components yields (m— 4)/2 independent eigenvectors of B[ with eigenvalues A^, and likewise for 
B' 2 . Hence the eigenvalues of B' are 

Afc = 2cos(— ), k = 3,5...m— 3, (A. 15) 

each with multiplicity 2, together with 0. If m = 2 mod 4, then is also in the set (|A.15p . and 
so the algebraic multiplicity of the zero eigenvalue is 3, even though its geometric multiplicity 
turns out to be only 2. However, this lack of full diagonalisability makes no difference to the 
computations of traces. 

We can now calculate lng B . Evaluating (|A.14|) using the eigenvalues just obtained, 

Det (/ - | (B') 2 ) sin 4 M s i n 4 M . . . sin 4 ^)ZL 



Det (I - \B 2 ) sin 2 in 2 3s:... gin 2 (rn^ 



m 

■ 2 3-n- ■ 2 5n -2 (m— 3)tt 

sin 2 ^ sin 2 *ZL • • • S in 2 ' 

mm m 



(A.16) 
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Labelling the final numerator X and the denominator Y, ()A.5|) at n = m and n = m/2 respec- 
tively implies that 

^ = -M^) 2 , y = ^- (A-i7) 



Hence 

X Detfl-i^') 2 ) 



Y ' Bet (I- \B 2 ) ~ ^msin 2 ^ 
and (|A.14|) agrees with (|4.16p . 



(A.18) 
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